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ABSTRACT 

We report the discovery of peculiar X-ray spectral variability in the binary radio millisecond pulsar 
PSR J0024-7204W in the globular cluster 47 Tucanae. The observed emission consists of a dominant 
non-thermal component, which is eclipsed for a portion of the orbit, and a thermal component, which 
appears to be persistent. We propose that the non-thermal X-rays originate in a relativistic intrabinary 
shock, formed due to interaction between the relativistic particle wind from the pulsar and matter 
from the main-sequence companion star, while the thermal photons are from the heated magnetic 
polar caps of the millisecond pulsar. At optical wavelengths, the emission exhibits large-amplitude 
variations at the orbital period, which can be attributed to heating of one side of the tidally-locked 
secondary star by the pulsar wind. The observed X-ray and optical properties of PSR J0024-7204W 
are remarkably similar to those of the low mass X-ray binary and X-ray millisecond pulsar SAX 
J1808. 4-3658 in quiescence. This supports the conjecture that the non-thermal X-ray emission and 
optical modulations seen in the SAX J1808. 4-3658 system in a quiescent state are due to interaction 
between the wind from a reactivated rotation-powered pulsar and matter from the companion star. 
The striking similarities between the two systems provide support for the long-sought connection 
between millisecond radio pulsars and accreting neutron star systems. 

Subject headings: pulsars: general — pulsars: individual (PSR J0024-7204W) — stars: neutron - 
X-rays: stars 



1. INTRODUCTION 

Radio millisecond pulsars (MSPs) form a separate class 
of rotation-powered pulsars, characterized by small spin 
periods, P < 25 ms, and spin-down rates, P — dP/dt ~ 
10 , implying surface magnetic field strengths of 

B oc (PP) 1 / 2 ~ 10 8 ~ 9 Gauss and characteristic ages of 
t = P/2P > 10 9 years. These neutron stars (NSs) are 
widely believed to be the end products of a more ex- 
otic channel of binary evolution, involving an extended 
period of accretion of matter and angular momentum 
from a close stellar companion in X-ray binary sys- 
tems l)Alpar et al.ll982UBhattacharva fc van den Heuvell 
11991(1 . After a period of 10 7 " 8 years, they acquire mil- 
lisecond spin periods and after accretion ceases are re- 
activated as radio pulsars. This theory is supported 
by the fact that over ^75% of the known MSPs in the 
Galactic disk are found in binary systems, whereas this 
is the case for only ~ 1% of the general pulsar popu- 
lation 1 . The strongest support so far for the connection 
between MSPs and X-ray binaries comes from the discov- 
ery of coherent millisecond X-ray pulsations in the low 
mass X-ray binary (LMXB) tran sient SAX J1808.4-3658 
l|Wiinands fc van der Klislll998D and sub sequently five 2 
other systems fsee iCampana et al"l ll998. for a review). 
However, a firm connection between the two classes of 
objects is yet to be established as no LMXB has been 
detected as a radio pulsar and no MSP has shown X-ray 
properties of a quiescent LMXB (qLMXB). 

1 Australia Telescope National Facility Pulsar Catalog, available 
at http:/ /www. atnf.csir.au/research/pulsar/psrcat 

s For the discovery of the latest X-ray millisecond pulsar see 
http: / /www.astronomerstelegram.org/?read=353 



The globular cluster 47 Tucanae (NGC 104, hence- 
forth 47 TucHs known to host at least 22 radio millisec- 

Jci 



ond pulsars liCamilo et alJl200(Tt IFreire et aUl2001l 1200.1 
Camilo & Rasiq l2005t and references therein). Many of 
these MSPs, however, have flux densities below the sen- 
sitivity of the Parkes radio telescope used to detect them 
and are only observable on rare occasions when strong in- 
terstellar scintillation focusing effects bring them above 
the detection limit. Thus, to date, 17 of the 22 MSPs 
have been detected sufficiently frequently to allow a pre- 
cise determination of their positions from radio timing so- 
lutions. During the only radio observation in which PSR 
J0024-7204W (henceforth 47 Tuc W or MSP W) was de- 
tected 3 , this 2.35-millisccond pulsar was found to be in a 
3.2 hour binary orbit and to undergo eclip ses but its po- 
sition could not be determined accurately l|Camilo et alJ 
2000). Fortunately, Hubble Space Telescope (HST) ob- 
servations of 47 Tuc have been used to match the binary 
period and phase of the modulations of a variable star 
with that of MSP W and thus determin e the position of 
the M SP with sub-arcsecond precision ((Edmonds et alJ 
2002) . Surprisingly, the optical counterpart of 47 Tuc W 
was found to be a >0.13 M Q main sequence star, which 
is unexpected for a MSP binary companion since in all 
but one other case the secondary star is greatly evolved. 
This intriguing result implies that either the companion 
is the star that very recently finished spinning up MSP 
W or that it is not the original secondary star but was in- 
stead exchanged for the low-mass (<0.1 M Q ) companion, 
which, in the process, was liberated from the binary. The 

3 After submission of this article for review, it was brought to 
our attention that this MSP has been detected again (P. Freire, 
private communication). 
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TABLE 1 

A SUMMARY OF THE X-RAY (Chandra) AND OPTICAL 
(HST) OBSERVATIONS USED IN THIS ANALYSIS. 



Telescope/ 


Epoch of 


Observation 


Instrument 


Observation 


ID 


Chandra/ACIS-S 


2002 Sep 29 


2735 


Chandra/ ACIS-S 


2002 Sep 30 


2736 


Chandra/ACIS-S 


2002 Oct 2 


2737 


Chandra/ACIS-S 


2002 Oct 11 


2738 


HST/ACS-HRC 


2002 Apr 5 


9028 a 


HST/ACS-HRC 


2002 Apr 12-14 


9019 


HST/ACS-HRC 


2002 Jul 24 


9443 


HST/ACS-WFC 


2002 Sep 30-Oct ll b 


9281 



a See also Edmonds et al. (2002). 
b Simultaneous with the Chandra 



ACIS-S observations. 



HST identification showed MSP W to be positionally co- 
incident with the X-ray source W29, first detected by a 
72-kilosecond Chandra X-ray Observatory ACIS-I obser- 
vation of the core of 47 Tuc lIGrindlav et alJl200l . The 
X-ray spectrum was found to be relatively hard, sugges- 
tive of non-thermal emission, atypical of the MSP sample 
in 47 Tuc, but simil ar to PSR J1740-5340 in the globular 
cluste r NGC 6397 (ID'Amico et al.ll200l IGrindlav et all 
2002), the only other known MSP with a non-degenerate 
stellar companion. 

2. OBSERVATIONS AND DATA ANALYSIS 

We have employed the superb spatial resolution of 
Chandra and HST to further investigate the unusual 
properties of the 47 Tuc W system. The X-ray dataset 
consists of four additional 65-kilosecond observations of 
the core of 47 Tuc performed with the Chandra ACIS- 
S instrument. These deep observations have resulted in 
the detection of ~300 sources with in the half-mass ra- 
dius of 47 Tuc ijHeinke et alJl2005^ including the X-ray 
counterparts for each o f the MSPs with known positions 
( Bogd anov et al]l2005lK The optical HST data consist of 
three sets of archival observations acquired with the High 
Resolution Camera (HRC) on the Advanced Camera for 
Surveys (ACS) and one set of observations performed 
contemporaneously with the Chandra observations us- 
ing the Wide Field Camera (WFC) on the ACS. Table 
1 summarizes all X-ray and optical observations used in 
this analysis. 

The initial data reduction and image processing of the 
X-ray observations were performed using the CIAO 4 3.0 
software package an d are discussed in detail elsewhere 
l)Heinke et alJ 1200*51) . For the spectral and variability 
analyses, the Chandra observations were first filtered to 
only include events with energies between 0.3 and 8.0 
keV. Subsequently, the X-ray photons from MSP W were 
obtained from a circular region with a 1" radius around 
the position of the optical counterpart as this circle en- 
closes >90% of the total energy for an assumed 0.5-1.0 
keV thermal source. The net count rate was then ob- 
tained by subtracting a background taken from three 
source-free regions on the image. To facilitate the spec- 
tral analysis, the ~300 net photons were grouped into en- 

4 Chandra Interactive Analysis of Observations, available at 
|http: / / asc.ha rvard.edu/ciao 



ergy bins, while requiring that each bin contain at least 
15 counts. For all spectral fits the hydrogen column den- 
sity to ward 47 Tuc was fixed at A^ H = (1-3 ± 0.3) x 10 20 
cm" 2 iGratton et alJl2003t iCardelli. Clavton. fc Mathisl 
H989UPredehl & SchmittJll995|) . An X-rav liehtcurve was 
obtained by folding the observations at the binary period 
using the latest ephemeris obtained from radio pulse tim- 
ing observations (Freire et al., unpublished data). 

The optical colors and variability of 47 Tuc W were 
first investigated in our ACS/ WFC data. Unfortunately, 
due to the large pixelscale of the WFC (0'.'049 per pixel) 
and the presence of several relatively bright stars around 
the position of 47 Tuc W, the faint optical source was not 
detected in the individual ACS observations. Since the 
images were dithered relative to each other with sub-pixel 
offsets, we could create combined images in each filter 
with improved spatial resolution using the Space Tele- 
scope Science Institute multidrizzle routines. In these 
phase-averaged images, the companion of 47 Tuc W is 
visible near the limit of detection. However, the large 
errors in the magnitude (^0.3-0.5 mag) prevented us 
from obtaining reliable lightcurves or a significant limit 
on the Ha emission expected from this system. Thus, 
we have used the higher resolution (0'.'027 per pixel) 
ACS/HRC archival data, where 47 Tuc W is clearly visi- 
ble when at maximum brightness. We measured mag- 
nitudes using aperture photometry on the distortion- 
corrected pipeline-processed images. A small aperture 
(2 pixel radius) was chosen to minimize contamination 
by the flux of nearby stars. For a few frames, mostly 
for those in which the counterpart of 47 Tuc W was 
near minimum brightness, a small change in the choice 
of background regions resulted in magnitude differences 
of up to 0.8 mag. We corrected these magnitudes to an 
aperture of 6 pixels using aperture photometry of sev- 
eral bright and relatively isolated stars within 200 pixels 
from the optical counterpart of 47 Tuc W. For each filter, 
we checked the variation of this correction as a function 
of time but found no indication for a systematic trend. 
This led us to apply a time-averaged aperture correction 
per filter. Next, w e used the filter-depe ndent encircled 
energy curves from ISirianni et all (|2005) 5 , that give the 
fraction of total source counts as a function of aperture 
radius, to co rrect the magnitude s to an infinite aper- 
ture (5'.'5, see ISirianni et al.ll2005l for details). We con- 
sider as unreliable those measurements for which these 
differences would result in flux differences of more than 
^30% and removed them from our dataset. The remain- 
ing data points do not show a systematic difference when 
compared to magnitudes derived for a different choice of 
background. Instrumental magnitudes were converted to 
magnitudes in the STmag system using the most recent 
filter-dependent sensitivities 6 . We used these fluxes to 
obtain sinusoidal fits to the lightcurves and make phase- 
dependent spectral energy distributions (SEDs). The 
times of mid-observation were converted to photomet- 
ric phases using the existing radio ephemeris (Freire et 
al., unpublished data). 



5 Available at http:/ /acs. pha.jhu.edu/instrument/calibration/photometry/ 

6 http://www.stsci.edu/hst/acs/analysis/zeropoints. 
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3. THE X-RAY SPECTRUM 

The X-ray spectrum of 47 Tuc W in the 0.3-8.0 keV 
range suggests that the emission consists of at least 
two components: a dominant non-thermal (power-law) 
and a fainter thermal component. The latter is fit 
equally well by either a simple blackbody (BB) or a 
NS hydrogen atmosphere (NSA) model ( Roman| 119871 
IZavlin. Pavlov, k, Shibanovll996ULlovdl2003|) . After ex- 
cluding the photons between orbital phases c/> = 0.4 — 0.7 
to correct for the bias in the best fit parameters intro- 
duced by the variability of the non-thermal component 
(see next section and Fig. 1) , we obtain the following best 
fit models. If we assume a BB model for the thermal com- 
ponent, we find that the spectrum can be described by a 
non-thermal component with a power law photon index 
of r = 1.13 ± 0.35 and a BB with an effective tempera- 
ture of T Gff = (1.56 ± 0.28) x 10 6 Kelvin and an effective 
emitting radius of R c s = 0.45 ± 0.43 km (xl = 1.22 for 
14 degrees of freedom), with a total unabsorbed flux of 
F x = (1.4 ± 0.3) x 10" 14 ergs cm~ 2 s" 1 . On the other 
hand, if we consider a NSA model and assume a gravita- 
tional redshift of z g = 0.31 at the NS surface, we obtain 
T = 1.14 ± 0.35, T eS = 0.95±g j4 x 10 6 K, R eS = - 2.3 
km (xl = 1.17 for 14 degrees of freedom), and a to- 
tal unabsorbed flux of F x = (1.5 ± 0.3) x 10~ 14 ergs 
cm~ 2 s _1 . In both cases, the non-thermal component 
contributes with ^75% of the total flux. F or a distance 
to 47 Tuc of 4.85 kpc l)Gratton et al.ll2003|) we obtain a 
total X-ray luminosity of Lx = (4 ± 1) x 10 31 ergs s _1 
and a non-thermal luminosity of Lx,nt ~ (2 — 3) x 10 31 
ergs s _1 (0.3-8.0 keV). All errors quoted above repre- 
sent ltr uncertainties in the derived parameters. We note 
that although an F-test does not indicate that a compos- 
ite model is statistically preferred over a pure power-law 
model (but with a steeper spectral index of T = 1.7±0.1), 
the X-ray lightcurve, discussed in §4, suggests that the 
thermal component in the spectrum is genuine. In addi- 
tion, the best fit values for T c g and R c g are quite sim- 
ilar to what is observed in MSPs with predominantly 
ther mal spectra, such as the majority of MSPs in 47 
Tuc l)Grindlav et alJEffll iBogdanov et al.ll2005fl . The 
marginally acceptable fits suggest that a second ther- 
mal component is probably present in the spectrum, as 
seen in the nearby MSPs J0437-4715 and JC I030+0451 
l|Zavlin et alJl2002t iBecker fc Achenbacljl2002l) . 

Thermal emission from MSPs is believed to origi- 
nate from the magnetic polar caps of the underlying 
NS, heated by a backflow of en ergetic particles from 
the pulsar magnetosphere (se e e.sr. lHarding fc Muslimovl 
120021 iZhang fc Chand I2003L and references therein). 
On the other hand, non-thermal emission from MSPs 
can be produced by two different mechanisms. In 
the most luminous MSPs, B1821-24 and B1937+21 
l)Becker k, Trumped Il999j) . the narrow X-ray pulse pro- 
files indicate that the non-thermal emission is beamed 
and must originate in the pulsar magnetosphere, near 
the light cylinder, and is probably synchrotron and cur- 
vature radiation from charged particles accelerated to 
relativistic speeds in the strong magnetic field. Alter- 
natively, non-thermal radiation can be produced by the 
interaction of the pulsar wind with the binary compan- 
ion or the interstellar medium. The resulting shock wave 
is expected to be a prominent source of high-energy 
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Orbital phase 

Fig. 1— (Upper panel) X-ray lightcurve of PSR J0024-7204W 
in the 0.3-6.0 keV band, folded over its binary period. Note the 
X-ray minimum at <j> fa 0.45 instead of <f> = 0.5, the position of the 
optical minimum. The uncertainty in the absolute phase is >10 — 4 . 
(Lower panel) The energy of each photon detected within 1" of the 
optical position of 47 Tuc W versus orbital phase. The data shown 
in this plot cover an integer number of binary orbits (21) to ensure 
uniform phase coverage. Note the remarkable absence of photons 
with energies above 2 keV and the reduction in the number of soft 
photons between orbital phases 0.4 and 0.7. Two orbital cycles are 
shown for clarity. 



photons, via the synchrotron or inverse Compton radi- 
ation produced by the highly accelerated particles at the 
shock front and in the flow downstream from the shock 
^Arons fc TWitfTflflfl. 

4. X-RAY VARIABILITY 

In the 47 Tuc W system, we do not know a priori the 
exact nature of the observed non-thermal emission. Al- 
though the Chandra ACIS-S observations preclude pul- 
sation analysis due to the limited time resolution (3.2 
s), the longer-timescale variability contains a wealth of 
information about the origin of the non-thermal X-rays. 
In particular, the lightcurve folded at the binary period 
exhibits dramatic variations as a function of orbital posi- 
tion (upper panel of Fig. 1). Based on Poisson statistics, 
we find that the observed variability is significant at the 
99.9% confidence level. 

Surprisingly, a comparison of the X-ray colors from 
different orbital phases shows that the spectrum of MSP 
W softens significantly around the minimum for ~30% 
of the orbit, with no counts detected above 2 keV (see 
lower panel of Fig. 1). In addition, at phases between 0.4 
and 0.7, we find a decline of about 70% in the number 
of soft photons (<2 keV). The remaining ~30% are con- 
sistent with being entirely due to the underlying thermal 
(BB or NSA) component although some of these pho- 
ton may be due to partial eclipses of the non-thermal 
source. In any case, there is strong indication that the 
non-thermal component is occulted for roughly one third 
of the orbit, most likely by the main-sequence secondary 
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star. Note that the underlying X-ray source cannot be 
eclipsed by diffuse material within or around the binary, 
as this would be manifested primarily by absorption of 
soft photons. 

The characteristics of the X-ray spectrum and 
lightcurve are suggestive of emission originating from two 
different locations within the binary. As stated above, 
the thermal component most likely originates from the 
heated polar caps of the MSP. On the other hand, the 
observed X-ray eclipses are inconsistent with a magne- 
tosphcric origin of the non-thermal photons. If these 
X-rays did originate from the pulsar magnetosphere, 
one would expect both emission components to undergo 
much shorter (at most 2-5% of the orbit) total eclipses, 
characterized by very rapid ingress and egress and cen- 
tered at exactly <j) = 0-5- The observed eclipses also 
suggest a localized non-thermal source, rather than an 
extended circum-binary X-ray emission region that could 
form, for instance, by interaction of the pulsar wind with 
the interstellar medium. 

The behavior of the hard X-ray emission can be nat- 
urally explained by the presence of a relativistic shock 
wave within the binary, produced by interaction between 
the pulsar wind and matter from the secondary star. 
This shock forms at some distance from the MSP where 
the ram pressure of the pulsar wind balances the pressure 
exerted by the gas from the companion star. The length 
of the X-ray eclipse (~30% of the orbit) suggests that 
this occurs much closer to the companion star than to the 
MSP. The shock can then be eclipsed by the secondary 
for a substantial portion of the orbit, even if the MSP is 
never occulted. As noted above, if the MSP is occulted 
by the secondary, which would occur for 80° < i < 90°, 
the eclipse should last for no more than 2-5% of the orbit. 
Unfortunately, the limited number of photons observed 
at the X-ray minimum (only 3 detected counts in a bin of 
width 10 minutes) prevent us from determining whether 
this is indeed the case and thus limiting the maximum 
allowed binary inclination. 

If the secondary star is underfilling its Roche lobe 
(RL), the observed shock could result from the pul- 
sar wind impinging on and evaporating material off 
of the surface of the companion or interacting with 
the stellar wind, generated due to the co-rotation 
forced by tidal coupling of the orbital motion to the 
deep convection layers o f the secondary star (see e.g. 
iWiiers fc Paczvriskilll993l and references therein). More 
likely, the main-sequence companion may be overflow- 
ing its RL, as a result of bloating due to irradia- 
tion by the energetic part i cle wind from the MSP 
(see e.g. IPodsiadlowskilll99lHD'Antona fc Erginl3ll993t 
iVilhu. Ergma. fc Fedorovalll994[) . The shock can then 
be produced by the interaction of the pulsar wind with 
a stream of gas from the companion issuing through the 
inner Lagrange point (Li). If this is indeed the case, 
it permits us to place a conservative lower limit on the 
orbital inclination by noting that at (f> ~ 0.45, the non- 
thermal X-ray source is completely occulted. For the RL 
geometry of a 1.4 M Q MSP and a 0.13-0.29 M compan- 
ion, correspon ding to the allowed r ange of inclinations 
90° > i > 27° ijEdmonds et al.ll200l . we find i > 35°. 

An interesting feature of the X-ray eclipse is the appar- 
ent sharp increase in the count rate immediately after the 
X-ray minimum. If real, this emission may be the signa- 




= 0.4 

Fig. 2. — A representation of the 47 Tuc W binary illustrating 
the key features of the system for the case of a Roche lobe filling 
companion, for assumed masses of 1.4 Mq and 0.15 Mq for the 
MSP and main sequence companion, respectively. Here, the orbit 
is depicted face-on (i = 0°) to better show the geometry of the 
X-ray emitting shocked region. The sense of rotation of the system 
is counterclockwise and is in the plane of the page. The dashed 
arrows point in the general direction of the observer roughly before 
(bottom) and after (top) the X-ray minimum. 



ture of Doppler boosting of the particle flow downstream 
from the shock as it t ravels around the secondary star 
l)Arons fc Tavani|ll993[) . This effect can also account for 
the possible drop in the countrate at <f> »s 1, correspond- 
ing to the portion of the orbit at which the particles 
in the shock are accelerated primarily in the direction 
away from the observer. Another notable property of the 
eclipse is its asymmetry, characterized by a rapid ingress 
and slower egress. In addition, the minimum in the X- 
ray flux occurs at <f> = 0.4 — 0.45, i.e. before the optical 
minimum at superior conjunction ((f) = 0.5). This sug- 
gests the presence of a swept-back shocked region that is 
elongated perpendicular to and is asymmetric about the 
semi-major axis of the binary. Such a geometry arises 
because as the gas from the secondary interacts with the 
pulsar wind it can no longer co-rotate along with the 
secondary star and must, therefore, fall behind. Figure 2 
illustrates the likely appearance of the non-thermal emis- 
sion region for the RL overflow case as inferred from the 
X-ray lightcurve. The geometry of the shock is similar 
if the companion star is contained within its RL except 
that for the case of evaporative mass loss the shocked 
"cometary" tail would originate at or very near the sur- 
face of the star. We note that a wind from the secondary 
star could also result in mass outflow through Li and 
a configuration similar to that shown in Figure 2. In 
all cases, when viewed at <f> ss 0.4, the apparent length 
of the shocked region is relatively small, explaining the 
short duration of ingress. Conversely, at (j> ~ 0.6 the pro- 
jected length of the emission region is larger, causing a 
gradual transition out of eclipse. The relative durations 
of ingress and egress as well as the length of the eclipse 
permit us to roughly constrain the geometry of the gas 
stream and the shock. Taking i — 90°, and a shock ap- 
proximately at Li, we find an upper limit on the length of 
the shocked stream of I < 2 x 10 10 cm, comparable to the 
radius of the companion, and a distance of d < 5 x 10 9 
cm, from the surface of the companion to the shock. 
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Fig. 3. — HST HRC light curves of the optical counterpart 
of 47 Tuc W in different filters folded at the photometric period. 
The flux is in units of 10 — 18 ergs s — 1 cm -2 A -1 . The data from 
observing programs 90f9, 9028, and 9443 are plotted with solid 
circles, gray circles, and open diamonds, respectively. For each 
filter the solid line is the sine curve with the best-fit amplitude 
and average flux level, with period and phasing from the most 
recent radio timing ephemeris (Freire et al., unpublished data). 
The dotted lines represent 2a error margins while the dashed lines 
give the average flux levels. The HST filter and A, the amplitude 
of variation (±lcr), are given in the upper left corner of each panel. 
Note the trend of decreasing amplitude with increasing wavelength. 
Two orbital cycles are shown for clarity. 



5. LARGE-AMPLITUDE OPTICAL VARIABILITY 

Figure 3 shows the HST ACS/HRC light curves in the 
F435W, F475W, F555W, and F606W bands. For these 
filters there is sufficient phase coverage to fit the am- 
plitude and average flux level of variation. For the fit 
we assumed a sinusoidal light curve and the latest radio 
ephemeris (Freire et al., unpublished data). To create 
SEDs, the fitted fluxes at maximum, mean, and mini- 
mum brightness were first corrected for the extinction 
towards 47 Tuc using the fi l ter-de pendent conversions in 
Table 32 of iSirianni et al J l)2005(l for the spectral type 
(G2) closest to that of the counterpart of MSP W (see 
below for the results of our blackbody fits). Since the 
detailed spectrum of the counterpart is unknown, it is 
difficult to estimate the corrresponding systematic er- 
rors. From comparison with the only two conversions 
provided for other spectral types (05 and MO), we ex- 
pect that these are small. The effective wavelengths that 
are used in the SEDs are derived with the task calcphot 
in IRAF/STSDAS 7 from folding the spectrum of a red- 
dened [E(B - V) = 0.023] 5500 K blackbody with the 
filter response curves. The resulting SEDs for photomet- 
ric phases 0, 0.25/0.75, and 0.5 were then fit to a BB 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 4. — The broadband spectrum of 47 Tuc W showing the 
X-ray data (circles) and optical data at maximum (squares) and 
minimum (diamonds). Also shown is the best fit model for the 
total emission (heavy black line), as well as for the individual emis- 
sion components: thermal emission from the MSP (dark grey line), 
blackbody emission from the secondary star at optical maximum 
and minimum (dashed and dot-dashed line, respectively) assuming 
Rett = 0.20 Rq, and non-thermal emission from the intrabinary 
shock (light grey line). The dotted lines delineate the bounds of 
the la- uncertainty in the best fit non-thermal (power-law) model. 



spectrum with fixed i? e ff but variable T c ff- 

We find that for the SED at maximum brightness, 
the data are best fit with T e g between —6800 K, for 
Rett = 0-15 R,0 (the radius appropriate for an unheated 
0.13 M star) and -5400 K for R cS = 0.26 R (corre- 
sponding to the RL radius for i — 90° and Mmsp = 1.4 
M Q ). Similarly, the fit to the SED at mean brightness 
results in a range of T ff ~ 6200 — 5000 K, while at min- 
imum brightness we obtain T e g ~ 5200 — 4400 K for the 
same range of stellar radii. We note that although the 
fits are consistent with R e s < Rrl, this does not neces- 
sarily mean that the companion is not RL filling. It is 
more likely that only a portion of the hemisphere facing 
the MSP is actually heated due to less efficient heating 
near the poles (see Fig. 2). Furthermore, small apparent 
values for R G ff can result from a viewing angle effect for 
i < 90°. Given the uncertainties in the inclination and 
the details of particle heating, we can only constrain the 
temperature difference between the heated and unheated 
areas of the secondary to AT C — 1000 — 2400 K, under 
the assumption that R c g at optical maximum is less than 
or equal to R e ff at minimum. 

An extrapolation of the non-thermal X-ray emission to 
optical wavelengths shows that the non-thermal emission 
produced in the intrabinary shock contributes negligibly 
to the optical emission, even if we take into account the 
large uncertainties in the absolute magnitudes of the op- 
tical data and the non-thermal X-ray photon index (see 
Fig. 4). Therefore, the observed large-amplitude vari- 
ability must be due to a temperature difference of the two 
hemispheres of the companion, which can be attributed 
to heating of one side of the tidally-locked secondary star 
by the MSP, and not from the shock. This scenario is 
consistent with the lightcurves in Figure 3, which show a 
trend of decreasing amplitude of the optical modulation 
with increasing wavelength. As we discuss later, both 
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the overall X-ray to optical spectrum and variability are 
very si milar to those observed in th e SAX J1808. 4-3658 
system l|Campana et alJl200ll200l . 

The combined X-ray and optical data permit us to 
determine the origin of the energy supplied to the sec- 
ondary star. In particular, we find that the total X-ray 
and UV flux from the intrabinary shock is insufficient 
to heat the face of the companion star to the observed 
temperature (T H ~ 5400 — 6800 K). Therefore, we con- 
clude that the heating of the secondary must instead be 
due to the relativistic particle wind of the MSP, gen- 
erated at the expense of the rotational energy of the 
underlying NS. The rate of this energy loss is given by 
E = 4tt 2 IP/P 3 , where / is the NS moment of inertia, 
typically assumed to be 10 45 g cm 2 . Unfortunately, 47 
Tuc W is rarely detected at radio frequencies which does 
not allow a determination of the pulsar spindown rate 
P. Hence, we have no direct measure of E. Nonetheless, 
we can use the correlation between the observed ther- 
mal X-ray luminosity (£x,t) and E found for the 47 Tuc 
MSPs, to obtain a crude estimate of this parameter. Us- 
ing Lx,t ~ (0.5 — 1) x I0 31 ergs s _1 and the empirical re- 
lation logL x ,T = 0.59 log £+10.0 l|Grindlav et all2002tt . 
we obtain £~lx I0 35 ergs s" 1 . Assuming an isotropic 
pulsar wind, we find that the total incident power on the 
secondary star is L m ~ 1 x 10 33 ergs s _1 . This value 
is sufficient to heat the face of the companion star to 
the observed luminosity (Lh ~ 1 x 10 32 ergs s -1 ), if we 
allow for a re-radiation efficiency of order 0.1. Further- 
more, the flux from the MSP wind at the surface of the 
secondary (~10 12 ergs cm -2 s _1 ) is conducive to bloat- 
in g of the star to a sign ificantly larger radius (see Fig. 1 
in lPodsiadlowskil[l99H) . provided that the particles con- 
stituting the wind deposit their energy below the stellar 
photosphere. Given that the pulsar wind is probably 
composed of highly relativistic electrons, positrons, and 
ions, which are able to penetrate deep into the stellar 
interior, this condition for expansion of the star can be 
satisfied, implying that the secondary star is most likely 
RL filling. Finally, the flux from the pulsar wind at Li is 
consistent with the observed value of Lx.nt ~ 3 x 10 31 
ergs s _1 , if we take int o account the particle accelera- 
tion efficiency e a ~ 0.2 l)Arons fc Tavanilll993l) . as well 
as the fraction of the wind interacting with matter from 
the companion (~ 10~ 3 ). 

6. EMISSION PROPERTIES OF THE INTRABINARY 
SHOCK 

The compact nature of the 47 Tuc W binary ensures 
that the intrabinary shock is located in a relatively strong 
pulsar magnetic field. Therefore, we expect synchrotron 
emission to be the principal energy loss mechanism in 
the shock. However, the proximity of the shock wave to 
the secondary star implies that the relativistic electrons 
in the shock are immersed in a "sea" of optical photons 
emanating from the secondary star. Thus, inverse Comp- 
ton scattering (ICS) may also be an important produc- 
tion mechanism for high-energy photons. 

The magnetic field strength immediately upstream of 
the shock is given by B 1 = [<r/(l + a^^iE/cfpd 2 ) 1 / 2 , 
where a is the ratio of the magnetic energy flux to the 
kinetic energy flux, f p is a geometric factor that defines 
the fraction of the sky into which the pulsar wind is 



emitted, and d is the distance between the MSP and 
the shock l)Arons fc Tav ani 1993). For simplicity, we will 
take d rj 6.4 x 10 10 cm, corresponding approximately to 
the distance from the MSP to Li assuming Mmsp = 1.4 
Mq and i = 60°. For an isotropic pulsar wind (f p = 1), 
we obtain B\ w 2 G and B\ w 30 G corresponding 
to the two possible cases of a kinetic energy dominated 
(a = 0.003 as in the case of the Crab nebula) and a mag- 
netically dominated (a 3> 1) wind, respectively. This 
implies a magnetic field strength beyond the shock of 
B 2 = 3i?i ~6G or B 2 ~ 90 G, respectively. 

Production of the observed e = 0.3 — 8 keV pho- 
tons via synchrotron, therefore, requires relativistic elec- 
trons with Lorentz factors 7 = 2.4 x 10 5 (e/i?2) 1 ^ 2 ~ 
1 x 10 4 — 5 x 10 5 . From these values, we obtain a ra- 
di ative loss time of £ SV nrh = 5 .1 x 10 8 (7-Bf) _1 ~ 1 — 60 
s ()Rvbicki fc Lightmanl 119791) . If we consider a power- 
law distribution of electron energies i.e. n e (-j) oc 7~ p , 
where y is related to the p hoton index by p = 2T + 1 
(Rvbicki & Lightman 1979), and assume a roughly cylin- 
drical shocked emission region with length I ~ 2 x 10 10 
cm and radius ^2 x 10 9 cm, we find that the required elec- 
tron density is n e ~ 10 7 or n e ~ 10 4 cm -3 , corresponding 
to a = 0.003 and a ^> 1, respectively. These values im- 
ply energy densities in electrons of U e ~ 10 2 — 5 x 10 3 
or U e ~ 10 5 — 5 x 10 6 ergs cm -3 . The value for 
a ^> 1 is comparable to the magnetic energy density 
U B = B%/8w ~ 1 - 300 ergs cm" 3 . Thus, based on 
equipartition of energy arguments, the wind is most likely 
magnetically dominated. 

We note that if ICS were to dominate the X-ray emis- 
sion, the radiating electrons need to be only weakly 
accelerated to 7 ~ 100, implying radiative loss times 
of hes = (tW^ph)isynch ~ 10 - 100 hours, where 
U p h = R c <tbTc I (i?*c) ~ 1 ergs cm -3 is the seed pho- 
ton energy density, while i? c , Tu , and i?* are the radius, 
temperature, and distance from the surface of the sec- 
ondary star, respectively. The implied U e ~ 10 7 ergs 
cm" 3 , however, greatly exceeds the values of Ub and 
[/ p h- Therefore, ICS is an unlikely soft X-ray production 
mechanism in 47 Tuc W. 

Future detailed multi-wavelength observations of the 
47 Tuc W system may reveal more information concern- 
ing the nature of the shock emission, which, in turn, may 
provide insight into the little understood properties of 
MSP winds, collisionless relativistic shocks, and particle 
acceleration mechanisms. The existence of a regularly 
eclipsed shock makes 47 Tuc W well suited for studies of 
these phenomena. 

7. DISCUSSION 

7.1. Comparison with Other MSPs 

The data presented here have revealed that 47 Tuc 
W exhibits X-ray variability that is unique among the 
MSPs (as evident in Fig. 1). These variations can be 
most easily explained by the presence of a relativistic 
shock within the binary that is regularly eclipsed by the 
secondary star. An intrabinary shock is also believed 
to be present in PSR B 1957+20, the so-called "black 
widow " eclipsing pulsar l|Fruchter. Stinebring. fc Tavlorl 
1988). This MSP has a non-thermal spectrum with 
T = 1.9 ± 0.5 and L x = 2.7 x 10 31 ergs s" 1 (0.5-7.0 
keV) l|Stappers et alJl20 031. for an assumed distance of 
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2.5 kpc l)Cordes Ez Lazioll2002j) . These values are compa- 
rable to those for MSP W and provide further support for 
the presence of a shock in the 47 Tuc W system. How- 
ever, we emphasize that 47 Tuc W is a fundamentally 
different binary than PSR B 1957+20 and similar eclips- 
ing systems, such as 47 Tuc J, O, and R (Bogdanov et 
al. 2005), which contain much less massive (~0.03 M Q ), 
greatly evolved secondary stars. 

MSP W posseses many characteristic similar to PSR 
J1740-5340 (J1740), the only known MSP in the globu- 
lar cluster NGC 6397 JP'Amico et al.l2001t iFerraro et all 
120031 ISabhi et al.M2003ft . Specifically. .71740 exhibits ra- 
dio eclipses, X-ray emission which is seem ingly hard and 
possibly variable ( Gri ndlav et alJ I2002D and is bound 
to a ~0.3 M Q "red straggler" companion that is cur- 
rently filling its Roche lobe. In addition, the shape 
of the Ha line observed in this 32.5-hour binary sys- 
tem indicates the presence of a swept-back stream of 
gas protruding from the companion, as the one pos- 
sibly present in 47 Tuc W (see Fig. 2). How- 
ever, unlike 47 Tuc W, J1740 shows no signature of 
wind-driven hea ting of the companion in the pho t omet- 
ric lightcurve llKaluzriv. Rucinski. fc Thompson! 120031: 
lOrosz fc van Kerkwiikll2003|) . This finding is consistent 
with the recently revised value of E « 3.3 x 10 34 ergs s _1 
l)Bassa fc Stappersll2004[) . which indicates that the flux 
from the pulsar wind incident on the ~1.6 R© compan- 
ion is too small to have a measurable effect. Therefore, 
in this system RL overflow is not the result of bloating 
by irradiation from the MSP but rather a consequence of 
evolution off of the main sequence. Perhaps the most pe- 
culiar feature of the J 1740 binary is the presence of He 
I absorption localized in a t hin longitudinal stri p near 
the surface of the secondary ijFerraro et alJ l2003). The 
shape of such a heated region could be due to irradiation 
of the companion by a highly anisotropic pulsar wind, 
preferentially emitted in the orbital plane. However, this 
would require the MSP spin and orbital angular momen- 
tum vectors to be exactly aligned. The thin strip can be 
more plausibly explained by the existence of an equato- 
rial wind, emanating from the tidally-locked companion 
star due to its forced rapid co-rotation, that is interacting 
with the pulsar wind. Alternatively, the He I absorption 
could originate in the swept back tail of material as the 
one shown in Figure 2 for 47 Tuc W. 

7.2. The MSP-LMXB Connection 

The unusual X-ray and optical properties of 47 Tuc W 
can serve towards finally establishing the long-suspected 
connection between LMXBs and MSPs. This is now 
possible because the 47 Tuc W system appears to be 
more typical of a qLMXB system, as it contains a main- 
sequence companion that may be RL filling, than a MSP 
binary system. In addition, the X-ray spectrum of 47 
Tuc W system exhibits remarkable similarities to that 
of the LMXB transient SAX J1808.4-3658 (henceforth 
J1808) in quiescence. XMM-Newton observations of the 
latter system have revealed that during the long pe- 
riods between outbursts (~2 years), this 2.01-hour bi- 
nary has an X-ray spectrum which is seemingly purely 
non-thermal, with T — 1.4i ' 3 and an X-ray luminos- 
ity of Lx = 5 x 10 31 ergs s _1 in the 0.5-10 keV band 
ijCampana et al.l|2002j) . Moreover, recent optical obser- 



TABLE 2 

A COMPARISON OF THE OBSERVED PARAMETERS FOR 47 
TUC W AND J1808. 



PSR J0024-7204W SAX J1808. 4-3658 



L x (ergs s 1 ) 


3 x 10 31 (0.3-8 keV) 5 x 10 31 (0.5-10 keV) 


r 


1.13 ± 0.35 (1<t) 


1.4±°j (3a) 


E (ergs s _1 ) 


~10 35 


? 


Li (ergs s -1 ) 


~1 x 10 31 


<3 x 10 31 


L irr (ergs s _1 ) 


~4 x 10 33 


>4 x 10 33 


Th (K) 


5400 - 6800 


6200t 8 3 ™ b 


L H (ergs s" 1 ) 


~1 x 10 32 


~4 x 10 31 


AT C (K) 


1000 - 2400 


~1000 



a L; is the estimated intrinsic bolometric luminosity of a 0.15 
Mq secondary star in the absence of heating from the pulsar; 
Th and Lh are the temperature and bolometric luminosity 
of the heated face of the companion. The other parameters 
listed are defined in the text. 

b S. Campana private communication. 



vations of the bloated ~0.05 M Q brown dwarf companion 
of J 1808 show sinusoid al variations at the orbital period 
ijCamnana et al.l2004D implying a temperature difference 
of AT C ^1000 ± 300 K between the two faces of the com- 
panion, quite similar to that reported here for 47 Tuc W. 
Finally, the observed non-thermal emission in J1808 also 
cannot account for the optical flux as well as irradiating 
luminosity [L- m > 4x 10 33 ergs s _1 ) required to produce 
the observed heating of the secondary star (see Fig. 1 of 
Cam pana et al.ll2004t) . The X-ray and optical properties 
of 47 Tuc W and J1808 are compared in Table 2. 

The primary difference between J 1808 in quiescence 
and 47 Tuc W is the lack of radio pulsations from the 
former. One possibility is that the radio beams are not 
favorably oriented, rendering J1808 undetectable at radio 
wavelengths. Alternatively, the radio signals from the 
nascent MSP expected in J1808 could be obscured by 
a diffuse circum-binary envelope of material, formed by 
the gas flowing out from the companion that is ultimately 
expelled from the system by the pulsar wind. Such an 
envelope is likely present around the 47 Tuc W system as 
well and may explain, in part, why this MSP is occulted 
for ~40% o f the orbit a nd is rarely detected at radio 
frequencies ijFreire 2005). If this is the case, it would be 
yet another common feature of the two systems. 

We propose that the great similarities in the values of 
the power law index and the non-thermal X-ray lumi- 
nosity derived for 47 Tuc W and J1808 point to a com- 
mon origin of the non-thermal X-ray emission in both 
systems. Also, it is highly probable that the source of 
energy supplied to the heated companion in both sys- 
tems is the same as well 8 . I f this is indeed the case, 
then, as original l y sugg ested bv lBurderi et alJ l)2003|) and 
iCampana et alJ (J2Q04), in the J1808 system there exists 
a pulsar wind, powered by the rotational energy of the 
NS, which is constraining the outflow of material from 
the companion and preventing its accretion onto the NS. 
As in the 47 Tuc W system, the presence of this relativis- 
tic particle wind is manifested via the non-thermal X-ray 

8 For J1808, unlike in the case of 47 Tuc W, no simultaneous 
X-ray and optical data are available and the best limit on Lx is 
<10 ergs s — 1 , as measured 2 days before the optical observations. 
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emission produced by interaction of this wind with mat- 
ter flowing out from the companion, as well as through 
the optical variability caused by heating of one side of 
the tidally-locked companio n. By analogy, th is would 
exclude an accretion disk (e.g lGarcia et alJ200lL and ref- 
erences therein) or a pro peller mechanism as the source 
of hard X-rays in J 1808 l)Wiinandsll2003t and references 
therein), since we know that in 47 Tuc W accretion is in- 
hibited by the pulsar wind. The remarkable similarities 
between 47 Tuc W and J1808 provide strong support for 
the claim that in episodes of quiescence, the NS in the 
qLMXB system J1808 has already been reactivated as a 
rotation-powered millisecond pulsar. This, in turn, gives 
indirect evidence in favor of the long-standing hypothesis 
that LMXBs are indeed progenitors of MSPs. 

7.3. X-ray Source Populations 

The implications of the 47 Tuc W system on studies 
of X-ray source populations are profound since low 
luminosity (Lx ~3x 10 31 ergs s _1 ) hard X-ray sources, 
which would otherwise be classified as cataclysmic 
variables (CVs), i.e. accreting white dwarf systems 
lIGrindlav et al.ll200l . may, instead, be MSP binaries in 
some cases. The same holds true for the low luminosity, 
hard spe ctrum q LMXBs categorized as black hol e (BH) 
systems llGarcia et al.ll200lt iRutledge et 101200011 based 
on the absence of thermal emission of a cooling NS. 
The optical heating effects of the pulsar wind, as seen 
in 47 Tuc W, would allow MSPs to be discovered 



and distinguished from CVs or BH qLMXBs since for 
both the low X-ray luminosities and the absence of an 
energetic wind do not result in significant heating of 
the secondary star for comparably short orbital periods. 
We note also that the discriminant of X-ray to optical 
flux ratios, Fx/Fy, which might be used to distinguish 
qLMXBs from CVs or MSPs is not useful for these 
transition qLMXB-MSP systems: J1808 and MSP W 
have log(Fx/-FV) ~ 0.7 and 0.4, respectively, which 
are typical values for CVs or qLMXBs but much larger 
than a "typical" MSP with an evolved companion. 
For comparison, the corresponding value for MSP U 
in 47 Tuc, whose companio n is likely a helium white 

dwarf llEdmonds et alJl200H) . is \og{F x /F v ) 1.15 

(IHeinke et al.ll2005l) . 
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